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We first show that the pions produced at high pr in heavy-ion collisions over a wide range 
of high energies exhibit a scaling behavior when the distributions are plotted in terms of a scaling 
variable. We then use the recombination model to calculate the scaling quark distribution just before 
hadronization. From the quark distribution it is then possible to calculate the proton distribution 
at high pr, also in the framework of the recombination model. The resultant p/n ratio exceeds one 
in the intermediate pr region where data exist, but the scaling result for the proton distribution is 
not reliable unless pr is high enough to be insensitive to the scale-breaking mass effects. 
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I. INTRODUCTION 

There are three separate and independent aspects 
about the hadrons produced at large transverse momen¬ 
tum ( pr) in heavy-ion collisions at high energies that col¬ 
lectively contribute to a coherent picture to be addressed 
in this paper. One is the existence of a scaling behavior 
at large pt that we have found by presenting the data 
in terms of a new variable. Another is the issue about 
the surprisingly large proton-to-pion ratio at moderate 
Pt (~ 2 - 3 GeV/c) discovered by PHENIX [I| in central 
AuAu reactions at ^fs = 130 and 200 GeV. The third 
issue concerns the hadronization process relevant for the 
formation of hadrons at large px and the applicability of 
the recombination model [ 3 ]. It is our goal to show that, 
in light of the scaling behavior of the 7r° produced, the 
recombination mechanism naturally gives rise to & p/n 
ratio that exceeds 1 in the 2 < px < 3 GeV/c range. 

Particle production in heavy-ion collisions at very high 
energies is usually described in terms of hydrodynamical 
flow Q, jet production at high px 3, thermal statis¬ 
tical model |5|, or a combination of various hadroniza¬ 
tion mechanisms 6j. In none of the conventional ap¬ 
proaches does one expect protons to be produced at 
nearly the same rate as the pions. If all hadrons with 
Pt > 2 GeV/c are regarded as products of jet fragmen¬ 
tation, then the known fragmentation functions of quark 
or gluon jets would suppress proton relative to pion by 
the sheer weight of the proton mass. Such a discrepancy 
from the observed data led some to regard the situation 
as an anomaly and proposed the gluonic baryon junction 
as a mechanism to enhance the proton production rate 
fj]. Their predictions remain to be checked by experi¬ 
ments. 

The parton fragmentation functions have been used 
even at low px in string models where the production of 
particles in hadronic collisions is treated as the fragmen¬ 
tation of diquarks, as done in the dual parton model )8|. 
There has been a long-standing dichotomy on whether 
particle production in the fragmentation region can bet¬ 
ter be described by fragmentation or recombination 


113 It is possible that the two pictures might be uni¬ 
fied in a more comprehensive treatment of hadronization 
in the future. Here we extend the recombination model 
to the central region at large px■ It should be recognized 
that an essential part of the recombination model is the 
determination of the distribution functions of the quarks 
and antiquarks that are to recombine. In the case of 
large-pT hadrons the underlying physics is undoubtedly 
hard collisions of partons and the associated radiation of 
gluons. If the parton distributions can be calculated just 
before hadronization, then the final step of recombination 
can readily be completed. If those distributions cannot 
be determined in pQCD, then the step between the ini¬ 
tiating large-pT parton and the resultant hadrons may 
efficiently be described by a fragmentation function, de¬ 
termined phenomenologically from experiments. Thus in 
that sense the two approaches, recombination and frag¬ 
mentation, are not contradictory, but complementary. 

We state from the outset that no attempt will be made 
here to perform a first-principle calculation of the parton 
distributions at large px before recombination. However, 
from the observed data on pion production in central 
AuAu collisions at the Relativistic Heavy-Ion Collider 
(RHIC), it is possible to work backwards in the recom¬ 
bination model to determine the quark (and antiquark) 
distribution at large px- On the basis of the quark distri¬ 
butions inferred, it is then possible to calculate the pro¬ 
ton distribution in the recombination model. The basic 
idea is that if there is a dense system of quarks and an¬ 
tiquarks produced in a heavy-ion collision whatever the 
dynamics responsible for them may have been (gluons 
having been converted to qq pairs before hadronization), 
then the formation of pions and protons (and whatever 
else) is prescribed by the recombination model without 
any arbitrariness in normalization and momentum de¬ 
pendence. 

One limitation of the recombination model as it stands 
at present is that it is formulated in a frame-independent 
way in terms of momentum fractions and is therefore in¬ 
applicable to a system where the particle momenta are 
low and the mass effects are large. The physics of recom- 
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bination is still valid at low momentum, but the details 
of the wave functions of the constituent quarks become 
important; they have not been built into the recombina¬ 
tion function that takes the simplest form in the infinite 
momentum frame. Thus our calculation of particles pro¬ 
duced at midrapidity is not reliable when pr is of the or¬ 
der of the masses of the hadrons under consideration. For 
protons we can trust the results only for pr > 3 GeV/c. 
For pions the lower limit of validity can be pushed much 
lower. 

Since our approach makes crucial use of the experimen¬ 
tal data on the pion spectrum as the input, it is essential 
to relate the spectra determined at different energies to 
an invariant distribution so that the scale-invariant re¬ 
combination model can be applied. To discover the ex¬ 
istence of an invariant distribution with no theoretical 
prejudices is a problem worthy in its own right. Fortu¬ 
nately, that turns out to be possible. The analysis for 
that part of the study will be presented below first to 
emphasize its independence from the theoretical model¬ 
ing of hadronization. It should be mentioned that the 
scaling of transverse mass spectra has been investigated 
recently 0 - The emphasis there has been on the depen¬ 
dences on the particle species and centrality for tut < 3.8 
GeV, while our focus is on the dependence on energy 
(17 < y/s < 200 GeV) for pr < 8 GeV/c. Thus the two 
studies are complementary to each other. 


II. A UNIVERSAL SCALING DISTRIBUTION 


The preliminary data of the pt distributions of n° pro¬ 
duced at RHIC at y/s = 130 and 200 GeV were shown 
by the PHENIX Collaboration at Quark Matter 2002 JIU 
for central AuAu collisions together with the WA98 data 
for PbPb collisions at y/s = 17 GeV 12]. They show that 
the level of the tail at large pr rises , as y/s is increased. 
We want to consider the possibility that the three sets of 
data points can be combined to form a universal curve. 

The 7r° inclusive distributions at midrapidity are in¬ 
tegrated over r] for a range of A rj = 1 so that the data 
points are given for the following quantity E2: 


f(PT,s) 


1 dN 

27 TPt dpx 


[ dr] ( 2np T N evt ) 1 d N ”° . ( 1 ) 
Jatj dprdr] 


In comparing the PHENIX data with those of WA98 one 
should recognize that in addition to the difference in the 
colliding nuclei there is a slight mismatch in centrality 
(top 10% for PHENIX and top 12.7% for WA98) Q- 
To unify the three data sets it is natural to first con¬ 
sider a momentum fraction variable similar to x f in lon¬ 
gitudinal momentum. However, so much momenta are 
taken by the other particles outside the Ar] = 1 range, 
it is unwise to also use y/s/2 as the scale to calculate 
the transverse momentum fraction. We assume that for 
every y/s there is a relevant scale K to describe the pr 
behavior relative to that scale. Let us define 


x = pr/K, 


( 2 ) 


and transform f(pr , s) to a new function $( 2 , K), where 
4yz,K) = K‘f(p T ,s) = /-A. (3) 

We adjust I\ for each s and check whether all three data 
sets coalesce into one universal dependence on z, which 
we would simply label as $(2), if it is possible. 

In Fig. 1 we show <f>(z), where the three symbols repre¬ 
sent the three data sets for the three energies. Evidently, 
the universality exists and is striking. While this behav¬ 
ior needs to be confirmed by more data, and the theo¬ 
retical implication remains to be explored, the existence 
of this scaling behavior is a significant phenomenologi¬ 
cal property of the pr distributions that suggests some 
underlying simplicity. It is like the KNO scaling of the 
multiplicity distributions P(n, s ) in pp collisions, where 
for y/s < 200 GeV they can be expressed by one universal 
scaling function with z = n/ (n) H 5 III 6 I]. 

The values of K that are used for the plot in Fig. 1 are 
in units of GeV: K = 1 (200), 0.9(130) and 0.717(17), 
the quantities in the parentheses being the values of y/s. 
The y/s dependence of K forms nearly a straight line, as 
shown in Fig. 2. Since the high and low energy data dif¬ 
fer both in colliding nuclei and in centrality, one does not 
expect strict regularity in how K depends on y/s. Nev¬ 
ertheless, an approximate linear dependence is a simple 
behavior expected on dimensional grounds. The straight 
line in Fig. 2 corresponds to the best fit 

K{s) = 0.69 + 1.55 x lO" 3 ^, (4) 

where y/s is in units of GeV. It should be recognized that 
the normalization of K(s) is arbitrary; it is chosen to be 
1 at y/s = 200 GeV for simplicity. If it is normalized 
to some other value at that point, the linear behavior in 
Fig. 2 is unchanged, only the scale of the vertical axis is 
shifted accordingly. The scaling property in Fig. 1 is also 
unchanged, the only modifications being the scales of the 
horizontal and vertical axes. Thus the absolute magni¬ 
tude of the dimensionless variable z has no significance. 

If the 2 dependence of <f>(z) in Fig. 1 were strictly lin¬ 
ear, so that it is a power-law dependence 

<F(z) oc z a , (5) 

then there would be no relevant scale in the problem. 
The fact that it is not a straight line implies that there 
is an intrinsic scale in the pt problem, which is hardly 
surprising. What is significant is that while there is no 
strict scaling in z, there is no explicit dependence on 
s. That is, at any energy we have the same universal 
function <I>(z), which will be referred to as the scaling 
behavior in s. That function can be parametrized by 

$( 2 ) = 1500 (2 2 + 2) ” 4 ' 9 , (6) 

which is represented by the smooth curve in Fig. 1. For 
large enough z Eq. © does have the form of the power 
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law given in Eq. 0 with a = 9.8. It is a succinct state¬ 
ment of the universal properties at high px- The depar¬ 
ture from Eq. © at small z reflects the physics at low 
Pt- Since there is no data on tt° for px < 1 GeV/c, 
the extrapolation of <f>( 2 ) to 2 < 1 is not reliable. How¬ 
ever, there is a more accurate determination of <f>(z) that 
includes the low z region when the charge 7r + data are 
considered; it is given in E2 , and is not needed here. 



FIG. 1: Scaled transverse momentum distribution of pro¬ 
duced 7r°. Data are from Ref. M H- The solid line is a 
fit of the data by Eq. 0 . 

Note that there is no fixed scale in px that sepa¬ 
rates the high- and low-px physics. Equation 0 gives a 
smooth transition from one to the other in the variable z, 
thus implying different ranges of values of the transition 
px at different s. 

While Eq. © gives a good parametrization of the scal¬ 
ing function $( 2 ) throughout the whole range of z, one 
notices, however, that the WA98 data at 17 GeV shows a 
slight departure from $( 2 ) at the high 2 end of that data 
set. It should be recognized that those data points have 
px > 3 GeV/c, which represents a huge fraction of the 
available energy at y/s = 17 GeV. In fact, one expects 
the violation of universality to be more severe at higher 2 
at that y/s, since energy conservation would suppress the 
inclusive cross section at higher px- What is amazing is 
that most of the WA98 data points are well described by 
$( 2 ), even though the corresponding px values take up 
a much larger fraction of the available energy than the 
other data points from RHIC. It demonstrates the signif¬ 
icance of the variable 2 in revealing the scaling property. 
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FIG. 2: The dependence of K{s) on y/s. The line is a linear 
fit. 


III. PION AND QUARK DISTRIBUTIONS IN 
THE RECOMBINATION MODEL 


Having found a scaling distribution for the produced 
7T° independent of s, we now consider the hadronization 
process in the recombination model in search for an origin 
of such a scaling behavior. In previous investigations the 
recombination model has been applied only to the frag¬ 
mentation region where the longitudinal momenta are 
large and the transverse momenta are either held fixed 
at low px or integrated over mu in We now consider 
the creation of pions in the central region of AA colli¬ 
sions and study the px dependence. Unlike the former 
case where the longitudinal momentum fractions of the 
partons are essentially known (from the structure func¬ 
tions), the px distributions of the partons in the latter 
case are essentially unknown. Indeed, it is the aim of this 
section to determine the parton px distributions from the 
7 r° distribution found in the previous section. 

Let us start by writing down the basic equation for 
recombination in the 3-space 


E- 


d 3 N n 


d 3 p 


d 3 pi d 3 p 2 
Ei E2 


F{pi,Pi)'R'n{pi,P2,V) 


( 7 ) 


where the left-hand side (LHS) is the inclusive distribu¬ 
tion of pion with energy-momentum ( E,p ). T(pi , p 2 ) is 






4 


the probability of having a quark at Pi and an antiquark 
at just before hadronization. TZ^ (p \, P2, p) is the in¬ 
variant distribution, Ed 3 N™/d 3 p 7 of producing a pion at 
given a g at pj and a q at p%. Note that lZ n has the 
dimension (momentum)^ 2 , same as the LHS. 

Writing the phase-space density in the form 


where 

F(zi,z 2 ) = K* f(p lT ,p 2T ) (14) 

R„(z 1 ,Z 2 ,z)=K { i z~ 2 . (15) 


d 3 p 

~E 


dy d(j>pr dpr, 


( 8 ) 


we define the inclusive distribution in pt, averaged over 
y and 0, 


d 3 N n 
Pt dp T 



( 9 ) 


where Ay is limited to one unit of rapidity in the central 
region. Our focus will be on the pr distribution at high 
Pt■ For the recombination distribution (p\ , p 2 , p) we 
need only consider the partons in the same transverse 
plane that contains p, since at high pt the partons with 
different yi are not likely to recombine. Indeed, we as¬ 
sume not only 2/1 = 2/2 = y, but also 0i = 02 = 0 
so that the partons and the pion are all colinear, and 
the kinematics can be reduced to that of a 1-dimensional 
problem. As in the usual parton model, the parton mo¬ 
mentum fractions in the hadron can vary between 0 and 
1, but the deviation in the momentum components of the 
partons transverse to the hadron p must be severely lim¬ 
ited because of the limited transverse size of the hadron. 
Thus we write 


Kn(puP2,P) = K -1 S (2/1 - 2/2) <5 ( 0 i - 02) 

8 ^ Vl + V2 -yj S 2 (p lT +p2 T - pt), (10) 

where 72.° is dimensionless, since 5 2 (pi T + P 2 T — Pt) has 
the dimension of lZ v (pi,p 2 ,p). If this 6 -function is fur¬ 
ther written in the colinear form due to the <5(0i — 0 2 ) 
in Eq. ([13 

$ 2 (Pi t +P2 t ~Pt) = s( ^ + ^ 2 - A — 

V 2 J p T 

Hpi t + P2 T - Pt) , (H) 


Since F(pi T ,p 2 T ) is the parton density in 
Pi T dpi T p 2 T dp 2 T , F(zi,Z 2 ) is the corresponding di¬ 
mensionless density in z\dz\Z 2 dz 2 ■ Equation m is 
now our basic formula for recombination in the scaled 
transverse-momentum variable. The total number of 
q and q is f dz\dz2 ZiZ2F(z±, Z2), which should be 
invariant under a change of scale 

z = Xx (16) 

so that 

x = pt/K', K’ = A I\. (17) 

The corresponding change on F(z\, Z 2 ) is that it becomes 

F , {x 1 px 2 ) = X 4 F(z 1 ,Z 2 ). (18) 

Thus the normalization of F{z\,Z 2 ) is scale dependent, 
as it should in view of Eq. Id- 

So far the recombination function R^izi, Z 2 , z) is not 
fully specified because 72.° has not been. In Eq. d 
the factor z 2 is associated with the dimension of the 
pion density, and the (5-function with momentum conser¬ 
vation. To introduce the pion wave function in terms of 
the constituent quarks, we rewrite Eq. d as 

R*(zi,Z2,z)=R°z- 2 G v (S 1 ,Z 2 ), (19) 

where i?° is a normalization constant to be determined 
and G-n-(£ 1 , £ 2 ) is the valon distribution of the pion 
Since the recombination of a q and q into a pion is the 
time-reversed process of displaying the pion structure, 
the dependence of J? ff (zi, Z 2 , z) on the pion structure is 
expected. During hadronization the initiating q and q 
dress themselves and become the valons of the produced 
hadron without significant change in their momenta. The 
variable in Eq. d denotes the momentum fraction of 
the ith valon, i.e., 


then Eq. 0 can be reduced to the ID form 


6 = Zi/z, 


( 20 ) 


dN v 

PTdpT 


J dpi T dp 2 T P\ T P 2 T F(p\tiPit)TQpt 2 


S 


f Pl T + P2t 

V PT 



( 12 ) 


where F(pi T ,p2 T ) is the qq distribution in pi T averaged 
over y and 0. 

We can reexpress this equation in terms of the scaling 
variable z = pt/K, introduced in Eq. |3- and obtain 


which is denoted by yi in the valon model mu , a no¬ 
tation that cannot be repeated here on account of the 
rapidity variables already used in Eq. d- In general, 
the valon distribution of a hadron h has a part specify¬ 
ing the wave-function squared, Gh , and a part specifying 
momentum conservation 

G h (Zi, ■■■) = Gfcfo, • ( 21 ) 


dN„ 

zdz 


= j dz 1 dz 2 z 1 z 2 F(zi,z 2 ) R- k (z 1 ,z 2 ,z) 


(13) 


where the functional form of Gh is determined phe¬ 
nomenologically. Although for proton G p is found to be 
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highly nontrivial |l9j . for pion G n turns out to be very 
simple 

G.(£i,6) = l, (22) 

which is a reflection of the fact that the pion mass is 
much lower than the constituent quark masses, so tight 
binding results in large uncertainty in the momentum 
fractions of the valons. Equation | 22 |) implies that the 
valon momenta of the pion is uniformly distributed in 
the range 0 < < 1 . 

What remains in Eq. m for us to determine is 
At this point we need to be more specific about the quark 
and antiquark that recombine. If the colors of q and q 
are considered, then the probability of forming a color 
singlet pion is 1/9 in 3 x 3. Similarly, for three quarks 
forming a proton the probability is 1/27 in 3 x 3 x 3. In 
the parton distributions, F qq for pion production involves 
two color triplets and F qqq for proton production involves 
three color triplets so the color factors work out just right 
in that the factors of 9 for qq and 27 for qqq are cancelled 
by the corresponding inverse factors in the recombina¬ 
tion probabilities. In other words, for the p/n ratio to 
be considered later, we can ignore the factors associated 
with the color degrees of freedom and proceed with the 
determination of F qq without specifying the quark colors 
and summing over them. 

The situation with flavor is not the same. For a uu 
pair and a dd pair, they can form 7 r° and 77 in the flavor 
octet. The branching ratio of 77 to 37 r° is 32.5% and to 
7 r + 7 T _ 7 r° is 22.6%. Thus for every 77 produced there is on 
average 1.27T 0 . Due to the higher mass of 77 we make the 
approximation that the rate of indirect production of 7r° 
via 77 is roughly the same as the direct production from 
uu and dd. If we now use qq to denote either uu or dd, but 
not both uu and dd, then each pair of qq leads to one 7r°. 
Since in a heavy-ion collision there are many quarks and 
antiquarks produced in the central region, it is reasonable 
to assume that the q distribution is independent of the q 
distribution so that we can write F qq in the factorizable 
form 

Fqq (zi,z 2 ) = F q (zi) F q {z 2 ), (23) 

where F q stands for either u or d distributions, and sim¬ 
ilarly for F q , but for ir° production q should be the anti¬ 
quark partner of q. 

The fact that we consider 77 production above, but not 
the vector meson p requires an explanation. We defer 
that discussion until the next section, after we have pre¬ 
sented the formalism for the production of protons. 

Returning now to the normalization of R^ (zi,z 2 ,z), 
we note that, using Eqs. frail, m and ( 12211 . 

J dzzR n {zi,z 2 ,z) = j ^ R°S ^ ~ 2 - lj 

= K (24) 

is the probability that a q at Z\ and a q at z 2 recombine 
to form a pion at any z. According to our counting in the 


second paragraph above, the total probability for qq —> 
7T° integrated over all momenta is 

J J fJ d zzR Aw) = l, (25) 

where Z is the maximum Zi, whatever it is. This normal¬ 
ization condition is scale invariant, and we find, using Eq. 
m ■ that 

K = 1- (26) 

Putting Eqs. ifTTitl - J23 and J 23 in 113 we obtain 

= J dz 1 dz 2 ^-^-F q (zi) F q (z 2 )8{zi + z 2 - z). (27) 

This is obtained from Eq. © where y and (f> are both ex¬ 
plicitly averaged over. The LHS is to be identified with 
$(z). Note that the 1/27T factors in Eqs. © and ©, 
where $(z) is defined, are there to render /(p T ,s) an 
average distribution in (/; that is the notation for the ex¬ 
perimental distribution, defined in G2- The distribution 
defined by us in Eq. © already includes the 1/27T fac¬ 
tor, so our dNn/zdz is just the experimental <f>(z). As 
we have mentioned earlier, the normalization of z has no 
significance. By means of a scale change in Eq. m we 
can move from z to x, or vice-versa, without changing 
the scale invariant form of Eq. (12711 . In Eq. © we found 
4>(z) to have the form 

$(z) = A(z 2 +c)~ n . (28) 

If we change z to x according to Eq. (trail . then by keeping 
the total number of pions invariant, i.e., 

dzzHz,K) = fdxz*'(x,K'), (29) 

we have 

$'0r, A K) = A 2 $(Ax, K). (30) 

It thus follows that 

$'(x) = A 2(1 - n) A (x 2 + c/A 2 ) ~ n . (31) 

Similarly, in the x variable the transformed quark distri¬ 
butions is 

F'(x 1 ,K') = X 2 F q (z 1 ,K). (32) 

Without having to specify the arbitrary scale factor A, 
let us work with the z variable and rewrite Eq. t?7i as 

< h(2) = J dzi Zi (1 - —) F q (zi) Fq(z - Zi). (33) 

We must now consider how the q and q distributions 
differ. Unlike the structure functions of the nucleon, 
where q and q have widely different distributions, we 
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are here dealing with the partons at high pt in heavy- 
ion collisions just before recombination. The dynamics 
underlying their pt dependences is complicated. Many 
subprocesses are involved, which include hard scatter¬ 
ing, gluon radiation, jet quenching, gluon conversion to 
quark pairs, thermalization, hydrodynamical expansion, 
to name a few familiar ones. At very large pt there 
are far more quark jets than antiquark jets, since the 
valence quarks have larger longitudinal momentum frac¬ 
tions than the sea quarks. By hard scattering the quarks 
therefore can acquire larger pt than the antiquarks. Thus 
in that way one would expect the pr distribution of the 
quarks to be very different from that of the antiquarks. 
However, that view does not apply to our problem. Those 
are the q and q that initiate jets, along with jets initiated 
by gluons. The conventional approach is to follow the jet 
production by jet fragmentation, which can be modified 
by the dense matter that the initiating partons traverse. 
As discussed earlier, our approach is not to delve into 
the dynamical origins of the q and q distributions, but to 
consider the recombination of q and q just at the point 
of hadronization. Such q and q are not the partons that 
initiate jets, but are the parton remnants after the hard 
partons radiate gluons which subsequently convert to qq 
pairs. Those parton remnants have similar momentum 
distribution for q and q , since gluon conversion creates 
q and q on equal basis; those partons are the ones that 
recombine to form hadrons. They are not to be confused 
with the jet-initiating hard partons that fragment into 
hadrons in the fragmentation model. In the recombi¬ 
nation picture those hard partons that acquire large pt 
immediately after hard scattering are not ready for re¬ 
combination; they lose momenta and virtuality through 
gluon radiation until a large body of low-virtuality quarks 
and antiquarks are assembled for recombination — a view 
that is complementary to the fragmentation picture. Of 
course, there are more quarks than antiquarks, since the 
number of valence quarks of the participating nucleons 
cannot diminish. For that reason we allow F q (z) to differ 
in normalization from F q (z). However, as a first approxi¬ 
mation we assume that their z dependences are the same. 

There is some indirect experimental evidence in sup¬ 
port of our assumption. In Ref. (jJ] the p/p ratio for 
central collisions is reported to be essentially constant 
within errors; more precisely, it ranges between 0.6 and 
0.8 for pt in the range 0.5 < pr < 3.8 GeV/c. Since 
p is formed by the recombination of three q , while p is 
formed from three q , a quick estimate of the q/q ratio is 
that it varies between 0.6 1 / 3 and 0.8 1 / 3 , i.e., from 0.843 
to 0.928. Such a narrow range of variation is sufficient for 
us to assume that F q (z) has the same z dependence as 
F q (z). For their relative normalization we take the mean 
p/p ratio to be 0.7. Thus we adopt the q/q ratio to be 

F~ q {z)/F q {z) = FL{x)/F' q {x) = 0.7 1 / 3 . (34) 

With this input we are finally ready to infer the quark 
distribution from the pion distribution. 
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FIG. 3: The solid line is the fit of the data as shown in Fig. 
1 (in a different scale), and the dashed line is the theoretical 
calculation of 4>(z) using the quark distribution in Fig. 4. 


We parameterize F q (z) by 

F q (z) = a (z 2 + z + z 0 ) m (35) 

and adjust the three parameters a, Zq and m to fit <E>(z) 
by using Eq. lEEl . We obtain an excellent fit with the 
values 

a = 90, zo = 1, to = 4.65. (36) 

In Fig. 3 we show in solid line the data represented by 
the formula in Eq. © and in dashed line the result of 
the theoretical calculation using Eqs. m~m- They 
coalesce nearly completely in the interval 1 < z < 8. The 
quark distribution F q (z) is shown in Fig. 4. To appreciate 
the pr range corresponding to z in Fig. 4, recall Eq. ©, 
Pt = zK , and Fig. 2 for K. Thus at yfs = 200 GeV, pt 
is z in GeV. Equations (1-351) and (® represent a main 
result of this study. What is important is that we have 
found a scaling quark distribution that is independent of 
s from SPS to RHIC, and perhaps to LHC. It is a succinct 
summary of the effects of all the dynamical subprocesses 
in heavy-ion collisions. The non-trivial z dependence in 
Eq. ©)j) indicates that there are intrinsic scales in the 
low-pr problem. 
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FIG. 4: Quark distribution in z. 


IV. THE p/iv RATIO 

The quark distribution obtained in the preceding sec¬ 
tion cannot be checked directly. Since it is the distribu¬ 
tion at the end of its evolution, massive dileptons would 
not be sensitive to it due to their production at the early 
stages. Proton production provides the most appropriate 
test, since liadronization occurs near the end. We shall 
therefore calculate the proton distribution at high pt and 
compare with the data on the p/it ratio. This is not a 
completely satisfactory venture, since the proton mass 
is large, so only at very high px can our scale invariant 
calculation be valid without explicit consideration of the 
mass effect. Present data on the p/n ratio do not extend 
beyond px ~ 3.8 GeV/c flj. Nevertheless, our calcula¬ 
tion should provide some sense on the magnitude of the 
rate of proton production at the high px end. 

The inclusive distributions in the scaled px variable 
can be obtained in the recombination model by general¬ 
izing Eq. tnu to the recombination of three quarks 

—— = [ dzidz2dz 3 Zi Z2 z 3 

zdz J 

F(zi,Z2,z 3 ) Rp(zi,Z2,z 3 ,z) (37) 

where F(zi, Z2,z 3 ) is given the factorizable form 

F(zi,Z 2 ,z 3 ) = F u (zi)F u {z 2 )F d (z 3 ). (38) 

As in Eq. 111)11 we relate the recombination function R p 


to the valon distribution, G p , of the proton 

R P {zi, z 2 ,z 3 , z) = R° z~ 2 G p (^ i,6,£ 3 ), (39) 

where G p has the general form given in Eq . <EU, and 
R p remains to be determined. In Ref. [l!j a detailed 
study of the proton structure functions has been carried 
out in deriving the valon distribution from the parton 
distributions that fit the deep inelastic scattering data. 
It is 


G p (ti,b,&)=g(Zi&) a ( 40 ) 

where 

a = 1.755, 0 = 1.05, (41) 

g = [B (a + 1, (3 + 1) B (a + 1, a, T 0 + 2)] . (42) 

Single-valon distributions G p (£i) can be obtained from 
the three-valon distribution by integration and are 
peaked around £ = 1/3, indicating that each of the three 
valons carries on average roughly 1/3 the momentum of 
the proton, their sum being strictly 1. Details of the 
valon model, described in [H|, are not needed for the fol¬ 
lowing. It is only necessary to recognize that the recom¬ 
bination of two u quarks with a d quark to form a proton 
has a probability proportional to the proton’s valon dis¬ 
tribution that accounts for the proton structure. The 
other point to bear in mind is that the valon distribution 
in the proton is obtained in the frame where the proton 
momentum is infinitely large so the finite masses of the 
proton and valons are unimportant. However, the valid¬ 
ity of that result when the proton momentum is only two 
or three times larger than its mass is questionable. With 
that caveat we proceed with our scale invariant calcula¬ 
tion and see what can emerge. 

As discussed in the preceding section, there is no need 
to consider the color factors for either pion or proton 
formation since hadrons are color singlets, but the flavor 
octets for these hadrons do introduce some factors. The 
\uud) state appears in 10 + 8 + 8' of 3 x 3 x 3; among 
them the first two contain A + and p. Thus the flavor 
parts of | (A+ \uud) | 2 and | (p \uud) | 2 are 1/3 for each. 
Since A + decays to p + 7r° and n + 7r + , | (jp | A + ) | 2 gives 
another factor 1/2. The spin decomposition of 2 x 2 x 2 
is 4 + 2 + 2, among which the A + component is 4/8 and 
p is 2/8. Putting the flavor and spin factors together, we 
have 


\(p\uud)\ 2 + |(p|A+) (A+ | uud) | 


_ 1 1 

“ 3 X 4 


1111 
- X - X - = -. 

3 2 2 6 


(43) 


We thus normalize R pi as we have done in Eq. (l25l) . by 


z 



dzzR p (zi,Z 2 ,z 3 ,z) 


1 

6 ' 


(44) 




In view of Eqs. m and (El we have 

r Z 3 


R, 


(ii A £A i 
L \ Z \z t z t z t J 6 


v L J-J- Z 

U 2=1 


(45) 


where z t = )T/ : Zi. Using Eq. El, the above integral can 
be transformed to 


„i 3 


m 


'U 


C t 


= 2.924 (46) 


with = Zi/Z and Q = JVCi- There is no explicit 
dependence on Z, and Eqs. ED and ED have been used 
in getting the numerical value in Eq. ED- It thus follows 
that 


R° p = 0.057. (47) 

At this point we should address the question why we 
consider A + production above, but not p production in 
the preceding section. For the production of 7 r°, if we 
are to consider the contribution from p ± (since p° does 
not decay strongly into 7r°), we would be extending our 
scope to other flavored states besides uu and dd. Then 
other vector mesons and higher resonances, such as K *, 
that can decay into 7r° must also be included. Similarly, 
for p production the consideration of other states beside 
uud would involve many resonances that can decay into p. 
The system is not closed without more phenomenological 
input beside 7r°. Thus for a closed system in which a 
prediction can be made, we limit ourselves to only the uu 
and dd in the meson states and uud in the baryon states; 
hence, only n°,ri,p and A + are considered. To include 
ud and du , we must also include uuu and udd , and so 
on. We surmise that if more resonances are included in 
both the meson and baryon sectors, the p/n ratio to be 
determined below would change somewhat; however, the 
result is not likely to differ by a factor greater than 2 . 

With the recombination function R p completely deter¬ 
mined, and the quark distribution F q ( Zi ) given by Eqs. 
(14511 and ED , we can now use Eq. ED to calculate the 
proton distribution in 2 . The result is shown by the solid 
line in Fig. 5, where only the portion 2 > 2 is exhibited. 
We have stated at the outset that the scale invariant form 
of dNp/zdz cannot be expected to be valid when the mass 
effect is important. The relevant value of z corresponding 
to the proton mass (let alone the A + mass) is 


l/2 

where itit = {m 2 + pip) and T p = 565 MeV. To con¬ 
vert this distribution to that for p we assume that only 
the normalization at pt = 0 needs to be adjusted. The 
p/p ratio at low pt is 0.6 pj. Since rriTdmT = Prdpr 
and the distribution in pt changes by a scale factor K 2 
given in Eq. EJ), where K = 0.9 for y/s = 130 GeV, 
the factor 4 in Eq. ED should therefore be changed to 
4 x 0.81/0.6 = 5.4. Expressing tut in terms of 2 by use 
of Eq. El with I\ = 0.9, we show the 2 dependence of 
the distribution for p in Fig. 5 by the short dashed line. 
The region 0.5 < 2 < 2 is left blank because our scaling 
result cannot be reliably extended into that region. Nev¬ 
ertheless, it is gratifying to observe that the theoretical 
calculation without any free parameters produces a pro¬ 
ton distribution at large 2 that is reasonable in normal¬ 
ization and shape and can smoothly be connected with 
the low ~pt distribution by interpolation. 



Z 


FIG. 5: Proton distribution in 2. Solid line is the theoretical 
result; the dashed line is the fit of data at low-py f20ll . 


z m = rn p /K, 


(48) 


which ranges from 1.3 at y/s = 17 GeV to 0.94 at 200 
GeV. As expected, the scaling violating effects are en¬ 
ergy dependent. Thus we should not regard the calcu¬ 
lated result to be reliable for 2 < 3. At very low pt the 
distributions of all hadrons can be given exponential fits 
in the transverse mass. The STAR data for most central 
collisions at y/s = 130 GeV J2(j give for p production for 
p T < 0.6 GeV 


1 d 2 N p 
rriT dvirdy 


4 exp [- (tot - m p ) /T p ] 


(49) 


With the proton distribution now at hand, we can cal¬ 
culate the p/tt ratio. For the pion distribution we use 
$(2) given in Eq. (0- For proton we use the calculated 
result based on Eq. ED- Their ratio, defined by 

Rph rW = ^/*(*) (50) 

is shown by the solid line in Fig. 6 . The preliminary data 
on the p/ 7 r ratio were reported in Ref. J, which we show 
also in Fig. 6 for both y/s = 130 and 200 GeV. Note that 
because it is a ratio there is no change in the normaliza¬ 
tions of R p / n for the two energies, but in transforming 
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from pt to z the factor I\ in Eq. 13 must be taken into 
account. Unlike the pion case the effects of the proton 
mass are not negligible for p T ~ 3 GeV/c, and one sees 
no scaling in s or z in Fig. 6. Our scale invariant calcu¬ 
lation is unreliable for 2 < 3 and shows a result that is 
obviously too high at 2 ~ 2. There seems to be a good 
chance that the theory and experiment can agree well for 
z > 4. In Fig. 6 we show two curves that can connect 
our scaling result with the data. The dotted curve is 
an eyeball fit of the 130 GeV data with a connection at 
^ = 3.5, while the dashed curve fits the 200 GeV data 
with a connection at z = 4. In the absence of a theo¬ 
retical study that takes the mass-dependent effects into 
account in the intermediate pt region, the only point we 
can make here is that it is not hard to produce a p/n 
ratio that exceeds 1 in the scale invariant calculation in 
the recombination model, but it does so in a region where 
both theory and experiment need refinement. Judging by 
what is self-evident in Fig. 6, we see no strong need for 
any exotic mechanism for proton production (as proposed 
in 0 ) beyond the conventional subprocess where three 
quarks recombine to form a proton. 



FIG. 6: Proton-to-pion ratio: solid line is the scaling distri¬ 
bution from calculation; data (preliminary) are from Ref. Q. 
The dotted and dashed lines are eyeball fits of the data as 
extrapolations from the scaling result. 


V. CONCLUSION 

The discovery of a scale invariant distribution $( 2 ) for 
pion production at intermediate and high values of pt 
in heavy-ion collisions ranging over energies in excess of 
an order of magnitude of variation is an important phe¬ 
nomenology observation that should be checked experi¬ 
mentally in great detail. Additional energy points should 
be added not only to strengthen the validity of the scaling 
behavior, but also to find the onset of scaling violation, 
if it exists. 

The phenomenological properties of hadron production 
provide useful insights into the hadronization process and 
into the nature of the quark system just before they turn 
into hadrons. The usual approach to the study of heavy- 
ion collisions is from inside out, following the evolution of 
the dense matter, either in terms of hydrodynamical flow 
or of hard parton scattering and subsequent hadroniza¬ 
tion by fragmentation [2l|. Our approach pursued here 
is from outside in, by starting from the observed scaling 
behavior of the pions produced and deriving the momen¬ 
tum distribution of the quarks that can give rise to such a 
behavior. That is accomplished by use of the recombina¬ 
tion model. There is no direct way to check the validity 
of the quark distribution thus obtained. However, we 
have used it to determine the proton distribution at high 
Pt where the mass effects are unimportant. The data 
on proton production have not yet reached that regime 
where the predicted scaling distribution can be checked. 
In the region where data exist on the p/7 r ratio we find 
that our calculated result, though not reliable, is in rough 
agreement with the imprecise data to the extent that the 
ratio exceeds 1, a feature that is notable. 

While the recombination model needs further work to 
take the proton mass into account at intermediate pr 1 its 
formulation in the invariant form has been developed here 
to treat the very high pr region. We have made the as¬ 
sumption that the quark and antiquark distributions are 
the same, apart from normalization, just before recombi¬ 
nation. That assumption is supported by the constancy 
of the p/p ratio in the PHENIX data in the central region. 
That experimental fact can also be used to lend credence 
to our general approach to hadronization that is treated 
as a recombination process, for which we have given ar¬ 
guments why the distributions of quarks and antiquarks 
should be similar before they recombine. In contrast, the 
fragmentation model would suggest a decreasing function 
of p/p in pt because of the dominance of quark jets over 
antiquark jets at large pt @1- 

In this paper we have only considered the energy de¬ 
pendence of the pt spectrum at fixed maximum central¬ 
ity. It is natural to ask what the dependence is on cen¬ 
trality. We have investigated that problem by making a 
phenomenological analysis of the data on centrality de¬ 
pendence without using any hadronization model, and 
found a scaling behavior very similar to what is reported 
here. The scaling distribution found there [TiJ includes 
the very small pt region in the fit, and is therefore more 
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accurate. But the fits in the intermediate and large pt 
region are the same. The implication on the centrality 
dependence of the p/ir ratio in the recombination model 
is still under study. 

To have an invariant quark distribution independent of 
s just before hadronization provides an unexpected pic¬ 
ture of the quark system. It suggests that the evolution 
of the system proceeds toward a universal form whatever 
the collision energy may be. We expect that universal 
form to depend on rapidity. The origin of such a scal¬ 
ing distribution in 2 is not known at this point and can 


form the focus of a program of future theoretical investi¬ 
gations. 
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